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(U)  The  radial  motion  of  gas  bubbles  in  water  driven  by  an  Intense  sound  field 
containing  two  closely  spaced  frequency  components  has  been  studied  via 
numerical  integration  of  the  equation  of  motion.  The  small  driving  amplitude 
perturbation  theory  solution  is  developed  for  comparison  with  the  numericeO. 
results.  Empheisls  is  placed  on  the  secondary  frequency  component  at  the 
difference  of  the  driving  frequencies  which  is  generated  by  the  nonlinear 
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Cases  where  the  natural  resonance  frequency  of  the  bubble  is  much 
greater  than,  approximately  equal  to,  and  less  than  the  driving 
frequencies  are  discussed.  Acoustic  driving  pressures  of  0.1  bar  to 
1.5  beir  were  used  in  eui  ambient  pressure  of  1.0  bar.  The  spectrum  of 
the  radial  motions  has  been  extensively  examined  and  many  examples  are 
shown  in  this  report. 

jill)  Agreement  between  perturbation  theory  predictions  and  the  numerical 

integration  is  good  for  the  first  and  second  order  terms  up  to  driving 
pressures  of  0.25  bar  in  all  cases.  All  cases  show  many  very  narrow 
lines  at  the  lower  driving  pressures . Broadband  noise  generated  by 
collapsing  bubbles  eventually  obscures  these  lines 

(u)  At  moderate  to  high  driving  pressures  the  difference  frequency  component 
is  larger  than  the  perturbation  theory  prediction,  especially  for  the 
case  of  the  driving  frequencies  being  below  resonance.  This  suggests 
that  larger  bubbles  driven  at  moderate  levels  may  be  the  most  useful 
in  enhancing  the  nonlinearity  of  water  although  tradeoffs  involving 
scattering  would  have  to  be  examined  in  each  particular  case. 
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I.  INTRODUCTION 


This  is  the  final  report,  under  Contract  N000l4-76-C-1095^  of  a 
study  of  the  response  of  small  air  bubbles  in  water  to  intense, 
mviltifrequency  acoustic  waves . The  approach  wa.s  to  numerically  integrate 
the  equation  of  motion  of  a bubble  using  a variety  of  bubble  sizes  and 

acoustic  driving  pressures . Because  the  nimierical  integration  reqid.red  , 

considerable  computer  time,  only  a limited  number  of  cases  was  studied. 

Basically,  this  is  a study  in  bubble  dynamics,  similar  to  studies  - 

1 2 ' 
conducted  by  Flynn  and  many  others . The  impetus  for  the  study,  however, 

came  from  the  field  of  nonlinear  acoustics.  Westervelt^  had  proposed  and  « 

Muir  and  Blue  had  confiimed  that  large  amplitude  sound  waves  containing  ^ 

two  frequencies  interact  in  water  to  produce  secondary  soxind  waves  at  the  ^ 

sum  and  difference  frequencies.  This  process  is  highly  inefficient.  ^ 

However,  if  microscopic  bubbles  are  introduced  into  the  water,  the  | 

5 8 

nonlinearity  of  the  medium  increases.  The  possibility  therefore  exists  j 

that  the  efficiency  of  nonlinear  generation  of  low  frequency  sound  can  be 
improved  without  seriously  degrading  all  of  the  process's  favorable 

features.^  While  experimental  work^'  has  verified  that  this  principle  , 

operates,  most  of  this  work  has  focused  on  acoustic  measxirements  and  little  ^ 

is  known  about  bubble  size  distributions  and/or  bubble  densities.  The 
present  study  was  undertaken  to  partially  bridge  the  gap  between  theory  and 
experiment . 

The  method  of  this  study  was  to  numerically  integrate  the  equation  of  ‘ 

motion  of  a spherical  bubble  in  a high  intensity  sound  field  consisting  of 

two  frequencies . Experimental  conditions  described  in  Refs . 10  and  12  , 

were  modeled.  While  the  motion  of  a single  bubble  will  not  yield  numbers 
which  can  be  directly  compared  to  the  experimenteil  data,  examination  of  the 
motions  of  bubbles  of  various  sizes  can  show  which  sizes  eu*e  important  in 

the  radiation  of  sound  at  the  difference  of  the  driving  frequencies.  i 1 
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This  study  determines  the  limits  where  perturbation  theory  is  valid 
euad  indicates  quantitatively  how  results  obtained  from  an  exact 
solution  differ  from  those  obtained  from  perturbation  theory  at  high 
pressures  in  complex  driving  fields. 

An  outline  of  the  theory  of  bubble  motion,  including  the  perturbation 
theory  solution  with  two  driving  frequencies,  is  given  in  section  II. 

The  numerical  integration  program  and  the  data  base  it  produced  are 
discussed  in  section  III.  The  results  axe  presented  in  section  IV, 
where  three  cases  are  examined  in  detail;  a summary  is  given  in 
section  V. 


II.  THEORETICAL  MCKGROIMD 


The  motion  of  a bubble  in  an  infinitesimal  amplitude  sound  field 
is  fairly  simple:  the  bubble  oscillates  at  the  driving  frequency. 

However,  the  case  of  bubbles  moving  under  a large  driving  pressure  is 
much  more  complex.  Above  a driving  pressure  threshold,  the  bubbles  can 
go  into  an  almost  complete  collapse  as  the  momentum  of  the  in- flowing 
water  overcomes  the  increasing  pressure  of  the  compressed  gas  within 
the  bubble.  During  these  collapses  the  radius  reduces  to  a very  small 
fraction  of  its  equilibrium  value . The  precise  form  this  collapse 
takes  depends  on  the  equation  of  state  of  the  very  hot  gas  within  the 
collapsed  bubble  and  is  not  well  understood. 

The  situation  with  a single  frequency  intense  sound  field  has  been 
studied  extensively  both  analytically^’ and  numerically The 
two- frequency  case,  however,  is  more  difficult  to  study  because  of  its 
complexity.  Analytical  methods,  other  than  in  the  small  amplitude  limit 
where  perturbation  theory  can  be  applied,  do  not  appear  very  fruitful. 
Therefore,  "the  differential  equation"  for  the  bubble  motion  has  been 
studied  numerically.  There  are,  in  fact,  several  differential  equations 
which  are  used.  The  four  most  common  ones^  are  shown  in  Fig.  1: 

(l)  Noltingk-Neppiras  - assumes  an  incompressible  fluid;  (2)  acoustic 
approximation  - based  on  allowing  for  a fixed,  finite  speed  of  sound; 

(5)  Herring- Flynn  - based  on  a better  accounting  of  energy  storage  on 
compression;  and  (4)  Kirkwood- Bet  he  - beised  on  specific  enthalpy  and 
tLLlowing  for  the  speed  of  sound  to  vary  with  position  and  time . These 
various  approaches  are  euLtemative  approximations  to  the  equation  of 
state  of  the  fluid.  All  give  essentially  the  same  result  except  in  the 
final  phase  of  bubble  collapse . To  study  the  broadband  high  frequency 
shock  waves  generated  by  bubble  collapse,  the  Kirkwood- Bethe  equation  is 
required.  In  addition  to  an  assimiption  about  the  equation  of  state  of  the 
fluid,  an  equation  of  state  for  the  gais  must  be  assumed.  Most  investigators 
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to  date  have  assumed  either  isothermal  or  adiabatic  conditions  for  the 
gas  with  a specified  ratio  of  specific  heats.  (Epstein  has  done 
some  work  with  adiabatic  cases  where  the  ratio  of  specific  heats,  7, 
takes  on  a second  value  during  the  finsil  stages  of  collapse . ) This 
study  hais  considered  only  adiabatic  motion. 

Whereas  this  study  concentrated  on  numerical  Integration  of  the 
acoustic  approximation  equation  and  not  on  analytical  methods,  a 
perturbation  solution  for  the  small  amplitude  case  is  useful  for 
comparison  purposes . The  differential  equation  considered  is 


where 

R is  the  bubble  radius, 
is  the  water  density, 

Pj^  is  the  pressure  at  the  bubble  wall, 
p„  is  the  ambient  pressure, 
c^  is  the  ambient  sound  speed, 
p is  the  viscosity  of  water, 

one  dot  indicates  differentiation  with  respect  to  time,  and  two  dots 
indicate  a double  differentiation  with  respect  to  time.  The  driving 
acoustic  field,  p , will  be  contained  in 

8i 

Pl  = Pb  Pa  ^ 

where 

p^  = p^  sinuo^t  + pg  sinu^t  , (5) 
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In  the  above  equations,  is  the  gas  pressure  in  the  bubble  and 
consists  of  the  compression  term  and  a term  involving  the  surface 
tension  o.  The  value  of  the  adiabatic  exponent,  y,  coxlLd  be  set 
arbitrarily  in  the  computer  program.  In  this  way  adiabatic 
(7=5/3)  or  isothermal  (7=1)  approximations  were  made.  is  the 
equilibrium  bubble  size  in  the  absence  of  an  acoxistic  sound  field 
and  the  driving  acoustic  pressures  at  angular  frequencies  and  cjd^ 
are  p^  and  p^,  respectively. 

To  obtain  the  simplest  infinitesimal  motion  approximation,  one 
can  examine 


I } 


and  let 


R 


-=l.x 

o 


(6) 


where  x is  a small  number.  If  only  one  driving  frequency  is  used,  then 
the  equation  in  x to  first  order  would  be 


X + 


57P„  sioDj^t 


X = 


(7) 


P R 
o o 


P R 
o o 


The  natural  bubble  resonance  frequency,  o)^,  in  this  approximation  is 


given  by 


37P„ 


° pr" 

o o 


(8) 


and 


X = Xq^  slnao^t  , 


(9) 


T7 


with 


'01  n 2.,  2 2v 

p R (OD  -COt  ) 
o o ' o 1 ' 


(10) 


For  two  driving  frequencies,  the  perturbation  theory  solution  to 
first  order  is 


X = Xq^  sincUj^t  + Xq2  sinojgt  , (ll) 

with  Xqj^  and  x^^  given  by  Eq.  (lO)  with  appropriate  subscripts.  This 
solution  may  then  be  used  in  an  approximation  to  the  equation  of 
motion  which  contains  terms  of  second  order. 


At  this  point  one  must  decide  on  the  relative  size  of  x^^  ^2  and 

R q/^2  2 ^2  2 wavelengths  at  2*  bubble  is  so 

small  that  is  larger  than  this  ratio,  the  terms  containing  R p /c 
01  o a'  o 

and  xR^c^  may  be  ignored;  otherwise  they  must  be  retained.  For 
simplicity  these  terms  will  be  dropped  since  we  are  looking  only  for 
an  illustrative  solution. 


Letting 


X = Xf  + Xs  , 


(12) 


x_  being  the  first  order  solution  and  x the  second  order  term,  the 

r S 

equation  of  motion  becomes 


(l+x)  X + OJ  X = 
' ' o 


^ (P^  sino^t  + P2  sinu^t) 


P R 
o o 


5 .2  2 (3.m)  2 

- 4-  X -in  • * X 
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The  term  in  braces  on  the  right-hand  side  of  Eq.  (13)  is  balanced  by 


x„  from  the  left-hand  side,  leaving 
F 


X + 03  X = 
S OS 


■*i^F  " 2 *F  % 


. (15) 


(U) 


or 


t 


1 


J t 
i ! 
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This  expression  contains  all  of  the  essential  features  of  a more 
complete  development  for  infinitesimeLL  waves . A simple  harmonic 
oscillation  is  driven  by  terms  at  0 Hz,  SUg,  co^-kOg,  auid 

This  perturbation  theory  development  is  incorrectly  handling  the  0 Hz 
term.  The  other  forcing  terms  should  be  meaningful  at  low  driving 

pressures  when  all  freqxiencies  are  far  from  resoneinces . The  terms 

*2  ••  / \ 

arising  from  the  R and  RR  in  Eq.  (1)  are  Included  here  but  were  absent 
from  the  analytical  development  by  RoUeigh  and  Clynch.^®  In  that 
development  the  bubbles  were  assumed  very  small  and  consequently  wsis 
assumed  much  larger  than  all  other  frequencies  in  the  problem.  In  the 
current  study  these  assumptions  were  not  made  and  therefore  these  terms 
have  been  kept. 


A diagram  of  the  frequencies  studied  is  shown  in  Fig.  2.  In  this 
diagram,  f^  is  the  bubble  resonance  frequency. 


3 

I 


T 


FIGURE  2 

FREQUENCIES  USED  IN  STUDY 

Difference  Frequency  f^  Center  Frequency  f^  i 

Bubble  Resonance  Frequency  f^  (Bubble  radius  In  brackets) 
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Not  n.l  1 combinations  of  difference  frequencies  and  center  frequencies  were 
processed.  The  largest  downshift  ratio  (f^/f^)  used  was  4.5,  with  a 
4.5  kHz  center  frequency  and  a 1 kHz  difference  frequency.  From  Fig.  2 
it  is  clear  that  the  approximation  of  large  is  valid  only  for  the 
20  pm  and  5 pm  bubbles . The  1000  pm  bubble  in  particular  does  not 
satisfy  this  criterion  since  it  was  actually  driven  below  resonance. 


In  all  cases  studied,  and  were  only  slightly  different  and  the 
driving  pressures  were  equal. . Therefore,  the  first  order  amplitude, 
and  will  be  approximately  equal  for  the  cases  studied  and  will  be 

much  greater  than  the  difference  frequency  second  order  amplitude. 
Examining  each  driving  frequency  in  Eq.  (l5)  separately,  one  finds  that 


(18) 


so 


(oj^+tOg) 


1 /PiPgX 


so 


1 /PlP2\  ^ - ^cu^csg 

nfl2\p  IT,  2 27,  2 2~  ^271" 

107  \ ^»/  (cs^  -o)^  )(a3^  -ojg  )(a)Q  -(cOi-a)g)  ) 


(19) 


(20) 


andXgo(2u>2)  can  be  obtained  from  Eq.  (l8)  by  substituting  ujg  for  o)^.  The 
equation  for  has  been  used  to  put  these  equations  in  the  form  of 
dimensionless  ratios.  While  these  equations  can  give  qualitative  results 
away  from  resonances  and  low  driving  pressures,  they  are  really  not 
applicable  at  the  acoustic  pressures  usually  present  in  experiments. 
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For  an  cu^  that  is  much  larger  than  aiiy  of  other  frequencies  involved, 

• these  equations  reduce  to  the  results  of  Rolleigh  and  Clynch.^*^  In  that 
limit  the  amplitudes  of  the  bubble  radial  oscillations  at  the  sum  and 
difference  frequencies  are  equal  and  are  tvice  the  value  at  each  second 
harmonic.  Thus,  in  the  leu-ge  limit,  the  amplitude  oscillations  at 
the  secondary  frequencies  are  partitioned  into  three  equal  parts:  at 
the  sum  frequency,  at  the  difference  frequency,  and  at  both  the  second 
harmonics . However,  this  partition  will  not  apply  to  the  radiated  sound 
field.  The  sound  radiated  will  be  proportional  to  the  second  time 
derivative  of  the  volume  oscillation,  causing  each  frequency  component 
of  the  sound  to  be  weighted  by  the  frequency  squared. 

In  the  general  case,  the  equations  are  quite  complex.  Even  though 
they  apply  only  at  small  acoustic  pressures,  they  can  serve  as  a guide  to 
frequency  dependencies.  In  addition  to  examining  absolute  levels,  the 

* ratios  of  x^^'s  to  the  values  computed  by  numerical  integration  will  be 
compared  to  provide  an  estimate  of  when  the  perturbation  theory  solution 
can  be  used  to  predict  spectral  shape. 
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III.  STUDY  PROCEDURE 


The  motion  of  a gas  bubble  in  a two- frequency  intense  soimd  field  was 
Investigated  through  numerical  integration  of  an  equation  of  motion. 

The  equation  integrated  by  the  computer  was  the  acoustic  approximation 
with  viscosity  added: 


The  output,  which  was  unevenly  spaced  in  time,  was  first  placed  on 
magnetic  tape  and  then  weis  interpolated  to  produce  eui  even  time  grid  which 
was  Fourier  transformed.  The  output  of  the  fast  Fourier  transforms  (FFT's) 
was  placed  on  magnetic  tapes  and  plotted.  (The  capability  to  plot  time 
histories  of  bubbles  also  exists.) 

A.  The  Integration  Program 

Equation  (2l)  was  integrated  by  program  GASBUBL5  on  the  CDC  5200 
computer  at  Applied  Reseeurch  Laboratories,  The  University  of  Texas  at 
Austin  (ARL:Ur).  This  program  accepts  control  cards  which  define  the 
characteristics  of  the  bubbles,  the  ambient  fluid,  the  driving  sound  field, 
the  initial  conditions,  and  the  time  interval  to  integrate.  The  equation 
of  motion  was  the  i Integrated  using  a fourth  order  Runge-Kutta  algorithm. 
This  integration  procedure  weis  chosen  because  it  allows  the  step  size  to 
vary.  The  outpu';  from  QASBUBL5  is  a magnetic  tape  containing  the  times, 
radii,  and  the  first  and  second  time  derivatives  of  R. 

The  integration  time  step  was  varied  automatically  in  the  program.  An 
absolute  maximum  time  step  of  one  hundredth  of  the  smallest  period  (bubble 
resonance  period  or  driving  sound  period)  was  used.  A convergence  criterion 
on  the  change  in  bubble  radius  was  applied  after  ecush  step.  If  this  test 


is  failed,  the  program  backed  up  one  step  and  halved  the  step.  If  the 
convergence  criterion  is  failed  twice  consecutively,  a major  restart  was 
done.  In  this  case  the  program  backed  up  20  steps,  reduced  the  time  step 
by  a factor  of  10,  and  continued.  A circular  buffer  of  the  100  previous 
points  was  kept  to  allow  for  five  overlapping  major  reductions.  Only 
after  a point  was  100  calculations  old  was  it  moved  from  the  circular 
buffer  to  storage  for  output  to  msignetic  tape.  Major  reductions 
in  step  size  were  also  triggered  by  a negative  bubble  radius. 

A problem  was  encountered  with  this  algorithm  in  the  case  of  5 Mm 

bubbles  undergoing  collapse . The  step  size  was  constantly  being  reduced 

to  the  point  where  accimiulated  round-off  error  made  the  calculation 

-90 

meaningless.  Therefore,  when  the  step  size  became  less  than  10  sec, 

-90 

the  step  size  was  reset  at  10  sec  and  kept  there . 

In  addition  to  being  reduced  automatically,  the  step  size  was 
automatically  doubled  if  20  steps  of  the  same  size  occurred.  If  the 
resulting  step  size  was  greater  than  the  absolute  maximum,  the  maximum 
was  used. 

1 15 

To  validate  the  program,  the  curves  given  by  Flynn  ’ were 
recel-culated . This  was  done  for  both  adiabatic  euid  isothemal  conditions . 
In  both  cases  the  GASBUBL5  curves  were  in  excellent  agreement  with  those 
of  Flynn. 

B.  Data  Base 

Center  driving  frequencies  and  difference  frequencies  which  matched 
experimental  data  on  hand  were  chosen  for  study.  A wide  range  of  bubble 
sizes  (5  Mm  to  1000  nm)  was  used  to  include  cases  where  the  resonance  was 
above,  nearly  equal  to,  and  below  the  same  driving  frequency.  Pressures 
of  0.1  bar,  where  the  perturbation  theory  should  hold,  to  1.5  bars,  where 
the  perturbation  theory  certainly  will  not  hold,  were  used. 


Ik 


Table  I Is  a summary  of  the  cases  that  were  integrated.  In  each 
case  marked  with  an  X,  the  full  set  of  pressures  was  run.  Each  integration 
with  two- frequency  data  was  run  out  to  a time  corresponding  to  100  difference 
frequency  periods . All  integrations  began  with  the  bubble  at  rest . 
an  cases  an  adiabatic  exponent  of  5/5  was  used. 


In 


SUMMARY  OF  DATA  BASE 


Items  marked  with  P are  partially  complete. 
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IV.  RESULTS  OF  THE  NUMERIC  INTEGRATIONS 

A thorough  analysis  of  all  the  Integration  results  has  not  yet  been 
made.  Three  Illustrative  cases  will  be  examined  which  shed  light  on  the 
motion  of  bubbles  in  a complex  high  pressure  sound  field.  In  all  cases 
an  adiabatic  equation  of  state  (7=5/3)  for  an  ideal  gas  was  used.  The 
ambient  pressure  was  taken  at  1.0  bar,  the  surface  tension  was 
75  dynes/cm,  and  the  viscosity  wajs  0.01  poise. 

All  the  cases  considered  had  a difference  frequency  of  400  Hz.  Two 
bubbles  were  used,  (l)  A 20  pm  bubble  with  a resonance  of  I85  kHz  is 
illustrative  of  the  small  bubble,  high,  frequency  resonance  case.  (2)  A 
500  pm  bubble  (7*13  kHz  resonance)  was  used  with  two  sets  of  driving 
frequencies:  one  where  the  frequencies  are  on  the  order  of  the 
resonance  frequency  (4.3  Wiz  and  4.7  kHz)  and  the  other  where  the 
frequencies  are  above  the  resonance  frequency  (15.8  kHz  and  l6.2  kHz). 

A.  Small  Bubbles  (f^»^,  f^) 

As  an  example  of  a small  bubble  case,  a 20  pm  bubble  driven  at 
4.3  l^z  and  4.7  kHz  is  considered.  In  the  present  context  "small"  is 
taken  to  mean  small  enough  that  the  resonance  frequency  (which  is 
inversely  proportional  to  the  equilibrium  radius)  is  much  greater  than  the 
other  frequencies.  In  this  case,  the  f^  is  I83  kHz  and  easily  satisfies 
this  criterion.  It  should  be  noted,  however,  that  the  surface  tension 
contribution  to  the  internal  pressure  at  equilibrium  is  only  0.075  bars. 
Therefore  this  is  still  a very  compliant  bubble. 

^ The  Fourier  transform  of  the  normalized  bubble  radial  motion,  x, 

i for  a drive  pressure  of  0.1  bar  per  primary  (5j^=5g=0.l)  is  shown  in 

f Fig.  3.  Here  5.=P./P  . The  most  obvious  feature  of  this  fig\ire  is 

^ 1 1'  00 

t ‘ the  neirrowness  of  the  lines.  In  fact,  the  lines  are  narrower  than  a 

I resolution  cell  and  a true  estimate  of  the  line  level  must  el-low  for 

the  fact  that  part  of  the  energy  shows  up  in  sidebsuids  for  any  finite  FFT 


i 
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if  the  line  is  not  exactly  in  the  cell  center.  A Hanning  window  has 
been  applied  to  these  transforms  to  minimize  this  effect,  but  it  is 
still  present. 


The  two  lines  of  the  drive  frequencies  appear  to  be  at  -37  dB 
from  a measurement  of  the  levels  in  Fig.  3-  A more  exact  value  using 
the  adjacent  cell  values  gives  -3^.4  dB,  in  excellent  agreement  with 
Eq.  (10)  which  predicts  -34.0  dB.  In  addition  to  the  lines  at  the 
driving  frequencies,  second  order  lines  at  the  difference  frequency  of 
400  Hz  and  at  frequencies  of  2fj^,  and  shown.  As 

predicted  by  Eqs.  (18)  and  (19)  and  discussed  in  section  II,  the  energy 
in  the  latter  three  lines  is  divided  evenly  between  the  two  second 
harmonics  on  the  one  hand  and  the  sum  frequency  on  the  other.  The 
spectrum  of  x on  an  expanded  freqxiency  scale  is  shown  in  Fig.  4.  Here 
third  order  mixing  effects  show  up  ais  lines  around  and  also  as 
sidebands  on  the  main  line  at  the  driving  frequencies.  The  exact 
general  perturbation  solution  for  these  lines  has  not  been  presented 
here,  but  in  the  limit  of  central  line  should  be 


60  log(l)  - 


= -85  dB 


which  appears  to  be  in  good  agreement  with  Fig.  4. 


A list  of  all  first  and  second  order  line  levels  and  the  perturbation 
theory  predictions  is  given  in  Table  II  for  these  driving  frequencies 
and  this  bubble.  Results  for  two  higher  driving  pressures  are  also  shown. 
The  perturbation  theory  resxilts  are  remarkably  accurate  at  5 of  0.25  but 
diverge  from  the  integrated  results  at  0.5  bar  driving  pressure. 

The  spectrum  of  the  0.25  bar  driving  pressure  is  shown  in  Figs.  5 
through  7.  The  higher  order  terms  are  present  through  many  orders  as 
shown  in  Fig.  6.  It  is  interesting  that  the  first  and  second  order  terms 
are  predicted  so  accurately  when  so  many  higher  orders  have  a significant 
amount  of  energy.  Figure  7 shows  this  spectrum  to  a much  higher  frequency 
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including  the  region  of  the  bubble  resonance.  In  this  case  a band  of 
energy  about  the  natural  resonance  frequency  exists,  but  no  narrov  line 
is  present.  The  first  peak  in  this  bend  is  at  106  kHz  while  the  last 
peak  at  the  same  power  level  is  at  223  kHz . This  broadband  feature  is 
caused  by  the  slow  variation  of  the  driving  frequency  pressure  being 
seen  as  a variation  in  the  ambient  pressure  by  the  20  pm  bubble  whose 
resonance  frequency  is  183  kHz  at  1 bar  ambient  pressure . This  effect 
will  be  discussed  later  in  section  IV. D. 

In  Figs . 8 through  10  the  spectrum  of  the  20  mn  bubble  driven  at 
0.5  bar  per  primary  is  shown.  The  number  of  lines  is  considerable,  but 
all  these  "low  frequency"  lines  are  still  sharp  in  Fig.  8.  Figures  9 
and  10  show  that  the  lines  form  "beat  groups",  separated  by  4.5  kHz,  which 
finally  merge  into  a general  background  level.  This  general,  background 
extends  to  well  beyond  the  natural  resonance  frequency.  Again  the 
natural  resonance  has  been  spread  by  an  effective  variation  in  ambient 
pressure . 

At  even  higher  driving  pressures  the  broadband  energy  from  bubble 
collapses  overwhelms  all  other  effects  and  even  the  driving  frequencies 
are  masked.  This  is  evident  in  Figs.  11  and  12,  which  show  spectra  of 
this  bubble  driven  with  1.0  bar  per  primary  signals.  In  Table  II,  the  lines 
at  fj^  and  fg  were  at  -I6.7  dB  at  the  0.5  bar  driving  level.  If 
Increased  by  only  6 dB  in  going  to  1.0  bar  (as  predicted  by  pertxirbation 
theory),  they  would  be  at  -10  dB,  which  is  15  hB  below  the  noise  level 
in  the  vicinity  of  f^  and  fg.  The  origin  of  the  broad  low  frequency 
peaks  around  550  Hz,  1100  Hz,  1200  Hz,  etc.,  is  unknown. 

B.  Medium  Size  Bubbles  f^<^+fg) 


To  illustrate  the  intermediate  case  a 500  pm  bubble  driven  with 
the  same  4.3  kHz  and  4.7  kHz  driving  frequencies  was  chosen.  The  natural 
resonance  frequency  of  the  500  pm  bubble  is  7*13  kHz. 


kHz 

Stic  Pressure 

0.50  bar 

0.50  bar 

3 

CO  o 

GO  U 

•—  c 

0) 

u 

c 

tc 

c 

o 

0) 

Q>  U 

•—  c 
^ 0) 

^ 3 N N 

1 J 

3 O’  Z Z 

CD  01  ^ ^ 

U 

ll-  CO 

u ^ ^ 

U) 

3 

o 

u 

cC 

CD 

3 M 
O’  U. 
O) 

u 


RADIAL  MOTION  SPECTRUM  OF  20  yiii  BUBBLE 


RADIAL  MOTION  SPECTRUM  OF  20  ym  BUBBLE 


The  spectrum  of  the  bubble  radial  motion  when  the  driving  pressure  is 
0.1  bar  per  primary  is  shown  in  Fig.  15.  In  this  case  all  lines  are  narrow, 
including  the  resonance  line  at  7. 15  kHz.  The  third  and  higher  order 
lines  are  also  quite  prominent  even  at  this  low  driving  pressure . Many 
of  the  frequency  components  on  Fig.  15  have  been  labeled.  The  order 
of  the  term  giving  rise  to  the  line  is  given  in  parentheses.  Some  , 

fourth  order  lines  show  up  on  this  graph. 

Spectra  for  this  case  at  driving  pressures  per  primary  of  0.25  bar, 

0.50  bar,  and  1.0  bar  are  shown  in  Figs.  14  through  16,  respectively.  At  ■ 

0.25  bar,  the  general  noise  level  is  increasing,  but  many  lines  can  still 

be  distinguished.  At  0.5  bar  and  1.0  bar,  all  lines  are  masked  by  the 

general  broadband  noise.  The  dip  formed  by  the  noise  decreasing  from  ’ 

0 Hz  and  the  increased  noise  from  the  broadband  (at  1.2  kHz  on  Fig.  l^) 

was  seen  in  many  cases  at  high  driving  levels.  An  interesting  feature  of 

Fig.  16  is  its  similarity  to  Fig.  11.  The  absolute  scales  sure  also 

approximately  the  same  when  the  normalization  by  is  removed.  This 

implies  that,  at  high  driving  pressures  in  the  two- frequency  case,  the 

bubble  spectrum  becomes  approximately  independent  of  the  original  bubble 

size.  The  period  between  successive  pressure  absolute  minimum,  which  is 

equal  to  the  difference  frequency  period,  apparently  determines  the 

maximum  size  to  which  the  bubble  can  grow.  If  this  expleuiation  is  correct, 

the  spectra  would  not  follow  this  pattern  when 

resonance  frequency . 

A summary  of  the  line  levels  measurable  in  the  0.1  bar  and  0.25  bar  j 

cases  Is  shown  in  Table  III.  As  with  the  small  bubble  case,  perturbation  ^ 

theory  predicted  the  0.1  bar  values  quite  well.  The  predictions  for  the  ! 

first  order  terms  also  held  true  at  0.25  bar  but  the  second  order  terms  i 

( 

are  in  error  by  up  to  5.7  dB. 
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C.  Large  Bubbles 


The  large  bubble  regime  was  investigated  using  the  same  500  urn 
bubble  described  in  section  IV. B.  with  driving  frequencies  of  15.8  kHz 
and  l6.2  kHz.  The  difference  frequency  is  therefore  the  same  as  in 
the  two  previous  sections. 

Figure  I7  shows  the  radial  motion  spectrum  of  the  500  ^lm  bubble. 

The  same  features  shown  for  the  medium  bubble  size  in  Fig.  15  are  apparent 
in  Fig.  17.  The  first  and  second  order  lines  are  narrow  at  the  levels 
predicted  by  perturbation  theory,  as  shown  in  Table  IV. 

The  spectra  of  the  same  bubble  driven  with  frequencies  of  15.8  kHz 
and  16.2  kHz  at  pressures  of  0.25  hai*>  0.50  bar,  1.0  bar,  and  1.5  bar 
are  shown  in  Figs.  I8  through  21,  respectively.  The  levels  of  the  lines 
can  be  followed  to  the  1.0  bar  driving  pressure  level  and  are  given  in 
Table  IV.  Only  at  the  1.5  bar  level  does  the  broadband  spectrum  emerge, 
and  even  then,  the  driving  frequencies  and  resoneuice  frequency  lines  are 
evident . 

The  perturbation  theory  values  Eire  in  good  agreement  through 
driving  pressures  of  0.5  bar/primary  with  the  exception  of  the  ^’^^''■^2 
which  is  1.5  dB  high  even  at  0.1  bar,  and  the  difference  frequency  line, 
which  is  beginning  to  grow  rapidly  at  this  level.  One  very  positive 
result  is  that  the  difference  frequency  line  is  much  larger  than 
perturbation  theory  predicts.  At  1.0  bar  it  is  about  20  dB  higher  than 
expected.  Although  its  precise  vadue  was  not  computed  by  the  program,  the 
400  Hz  line  is  still  evident  above  the  noise  in  Fig.  21  and  probably  has 
a value  around  -23  dB. 

D.  Effective  Resonance  Freqiiency 

In  several  of  the  previous  figures  (Fig.  7 for  example)  the  resonance 
WEIS  spresul  out.  This  occurred  when  the  resonEuice  frequency  weis  higher 


Perturbation  Theory 


{ 
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than  the  driving  frequency  and  the  total  driving  pressure  approached  or 
exceeded  the  ambient  pressure.  An  example  of  this  type  of  spectrum  is 
shown  in  Fig.  22.  In  this  case  the  20  pm  bubble  had  been  driven  with 
only  one  frequency  to  make  the  interpretation  clearer. 

The  natural  resonance  frequency  of  the  bubble  is  I83  kHz  at  a 
static  pressure  of  1 bar.  The  driving  pressure  was  1 bar  at  2 kHz. 

During  the  maximum  pressure  the  resonance  would  therefore  be  X I85  kHz 
or  259  During  the  negative  part  of  the  acoustic  pressure  the  total 

pressure  will  be  zero  and  the  resonance  frequency  will  be  zero.  In 
Fig.  22  one  sees  that  the  spectrum  starts  at  0 Hz  but  extends  to  about 
290  kHz  or  1.6  times  the  ambient  pressure  value. 

The  time  waveform  of  the  bubble  motion  is  shown  in  Fig.  23  along 
with  the  acoustic  driving  pressure . One  sees  that  the  bubble  expands 
to  many  times  its  equilibrium  radius  during  the  negative  pressure  phase, 
acquiring  energy  which  will  show  up  at  the  bubble  resonance  frequency 
once  the  bubble  begins  to  oscillate.  Because  the  bubble  frequency,  even 
on  the  first  few  cycles,  is  much  higher  than  2 kHz,  the  bubble  experiences 
a constantly  increasing  ambient  pressure  as  it  goes  through  each 
oscillation  and  increases  its  frequency  until  the  pressure  begins  to 
decrease.  At  that  time  the  oscillations  ajre  subsumed  in  the  general  bubble 
oscillation  expansion  and  the  cycle  repeats.  The  bubble  is,  in  effect, 
executing  an  IM  slide. 

It  should  be  noted  that  the  motion  is  cyclic  even  though  there  is 
no  rational  fraction  relating  the  "natural"  resonance  at  I83.2  kHz 
and  the  driving  frequency  (2  kHz). 
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SUMMARY 


The  equation  of  motion  of  a gas  bubble  in  water  has  been 
numerically  integrated  for  the  case  of  being  driven  by  two  intense 
acoustic  waves  at  frequencies  f^  and  fg.  Emphasis  has  been  on  the 
secondary  frequency  components  generated  by  the  nonlinear  motion  of 
the  bubble,  particularly  the  difference  frequency  cases 

have  been  examined  in  detail:  driving  frequencies  far  below  resonance 
(small  bubble),  driving  frequencies  neax  resonance,  and  driving 
frequencies  above  resonance  (large  bubble). 

The  numerical  results  were  compared  to  results  obtained  by 
perturbation  theory.  Agreement  was  excellent  at  driving  pressures  up 
to  one  quarter  of  the  ambient  pressure  in  the  small  bubble  case,  while 
the  near- resonance  results  had  already  begun  to  diverge  at  this  level. 

The  large  bubble  results  were  in  good  agreement  with  perturbation  theory 
at  half-ambient  driving  pressures  except  for  the  difference  frequency 
component,  which  was  3.5  dB  larger  than  predicted. 

In  all  cases,  at  medium  pressures  a myriad  of  sharp  lines  were 
present,  sometimes  resulting  from  fourth  order  interactions.  At  high 
driving  pressures  (equal  to  or  greater  than  ambient)  sharp  bubble  collapses 
created  broadbeuid  noise  which  obliterated  the  sharp  lines . This  indicates 
that  driving  pressures  below  the  cavitation  level  may  be  more  useful  in 
generating  difference  frequency  sound. 

The  rather  large  enhancement  of  difference  frequency  level  over  the 
perturbation  theory  solution  in  the  large  bubble  case  indicates  that 
large  bubbles  are  probably  more  useful  than  small  ones . In  addition  the 
lines  are  present  to  much  higher  driving  pressure  levels.  In  designing 
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experiments  and  equipment  the  problem  of  scattering,  which  increases 
as  the  fourth  power  of  the  radius,  must  be  considered  and  a detailed 
trade-off  study  would  be  required. 
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APPENDIX 

MAGNETIC  TAPE  FORMATS 


A.  GASBUBL5  Output 


GASBUBL5  outputs  four  multiplexed  clmnnels  of  data  to  magnetic  tape 
in  SPIO  format.  These  cheuinels  are  time  (T),  bubble  radius  (R),  bubble 
wall  velocity  (R),  and  bubble  wall  acceleration  (K)  in  that  order.  The 
data  are  grouped  in  units  of  2048  points  per  array.  Each  of  these  units 
is  called  a "sequence"  and  consists  of  two  physical  records:  an  identi- 
fication record  (IDR)  eind  a data  record  (DR).  The  IDR's  are  121  CDC  3200 
integer  words  long  and  the  DR's  are  8096  (4*2048)  CDC  5200  Integer  words 
long.  The  records  were  written  with  BUFFER  OUT  in  odd  parity  at  80O  bpi. 
Each  sequence  takes  3-6  ft  of  tape  allowing  about  640  sequences  to  be 
placed  on  a 2400  ft  tape. 

Each  IDR  contains  coded  information  to  allow  easy  l/O  using  routines 
available  at  ARL:Ur  and  information  about  the  parameters  used  in  the 
integration.  Each  IDR  is  equlvalenced  in  GASBUBL5  to  an  integer  (ID) 
and  a floating  point  (XD)  array.  Because  it  requires  two  24-bit  CDC  5200 
words  to  store  one  48-bit  floating  point  word,  XD  is  only  60  elements 
long.  The  storage  format  of  the  IDR  is  given  in  Table  AI. 

The  four  multiplexed  data  channels  have  been  scaled  and  stored  as 
24-blt  integers  prior  to  being  placed  on  meignetic  tape.  Tlie  scale 
factors  and  other  quantities  necessary  to  retrieve  the  true  values 
of  the  data  are  stored  in  the  IDR.  Because  of  automatic  scaling  within 
GASBUBL3,  these  scale  factors  may  not  be  the  same  from  sequence  to  sequence. 
To  recover  the  data  on  a CDC  3200  computer,  one  must  input  the  IDR  and  set 


RZERO 

= XD(52), 

SCLT 

= XD(42), 

SLCR 

= XD(43), 

SCLRD 

= XD(44),  and 

SCLRDD 

= XD(45). 

*A  detailed  explanation  of  this  format  ia  given  in  Ref.  I8.  However,  the 
information  in  this  appendix  should  be  sufficient  for  a programmer  to  be 
able  to  read  the  tapes. 


TABLE  AI 

GASBUBL5  IDENTIFICATION  RECORD 
Integer  Array 

ID  Index  Storage  Usage 

1 Number  of  elements  in  each  array  In  following  DR  (2048). 

2 Sequence  number,  a seqi;iential  record  index. 

11  Initial  suggested  sample  frequency  (kHz) 

19,21  Used  by  l/O  routines. 

91  Number  of  steps  to  back  up  in  major  restart  (20). 

93  Nimiber  of  this  sequence  within  current  integration.  If  negative, 

indicates  last  sequence  in  integration. 

9^  Maximum  number  of  consecutive  step  size  reductions  before  doing 
a major  restart  (l). 

100  Number  of  aj^rays  multiplexed  in  DR  (if). 

121  Characters  "SPIT",  used  to  locate  ID  records. 


Floating  Point  Array 


XD  Index 

Storage  Usage 

22 

7 

Adiabatic  exponent 

25 

a 

Surface  tension  (dynes/cm) 

2k 

P 

0 

Ambient  pressure  (bars) 

25 

M. 

Viscosity  (poise) 

26 

Po 

Fluid  density  (g/cm^) 

27 

^1 

First  driving  frequency  (kHz) 

28 

Pi 

Amplitude  of  pressure  at  frequency  f^^  (bars) 

29 

^2 

Second  driving  frequency  (kHz) 

30 

Pg 

Amplitude  of  pressure  at  frequency  fg  (bars) 

31 

9 

Beginning  phase  difference  in  driving  amplitudes  1 

32 

R 

0 

Equilibrium  radius  (pm) 

33 

P 

V 

Vapor  pressure  inside  bubble  (bars) 

34 

t 

Initial  time  for  this  integration  (psec) 

35 

R 

Initial  bubble  radius  for  this  integration  (pm) 

36 

R 

Initial  bubble  wall  velocity  for  this  Integration 
(pm/psec=m/sec) 

Numbers  in  parentheses  axe  default  values . 


5*^ 


T 
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XD  Index 


TABLE  AI  (Cont’d) 


Fractional  part  of  period  which  step  si7.e  will  not  exceed  (O.Ol) 

Factor  by  which  time  step  increased  and  reduced  in  minor 
adjustments  (2) 

Sj^  Test  factor  on  Runge-Kutta  k's  (O.Oj) 

f^  Bubble  resoneince  frequency  (kHz) 

C Sound  speed,  if  zero,  acoustic  terms  omitted  from 
® integration 

SCT  Scale  factor  for  time  array 

SCR  Scale  factor  for  R axray  (this  DR) 

SCR  Scale  factor  for  R array  (this  DR) 

SCR  Scale  factor  for  R array  (this  DR) 

T(1)  Beginning  time  in  this  DR  (usee) 

R(1)  Beginning  radius  in  this  DR  (pm) 

R(l)  Beginning  bubble  wall  velocity  in  this  DR  (^im/^lsec) 

T(liast)  Last  time  in  this  DR  (usee) 

R(Last)  Last  radius  in  this  DR  (pm) 

R(Last)  Last  bubble  wall  velocity  in  this  DR  (pm/usec) 

Atj Maximum  allowed  At  (usee) 

At^ln  Minimum  At,  this  DR  (^xsec) 


At  Maximvmi  At,  this  DR  (^sec) 
max  ’ ' 

Program  name  (Characters  "GASBUBL5") 


I 

I 


The  DR  can  then  be  read  into  an  integer  array  (IA)  dimensioned  (h,20k6). 
The  time  in  microseconds  is  given  by 

T(I)  = IA(1,I)  * SCLT  , 
the  radius  in  micrometers  by 

R(l)  = IA(2,I)  * SCLR  * RZERO  , 

R in  micrometers  per  microseconds  by, 

RD(I)  = IA(3,I)  * SCLRD  * RZERO  , 

and  R by 

RDU{I)  = IA(4,l)  * SCLRDD  * RZERO 
B.  Resampled  Output 

The  program  GASRSMPL  accepts  the  output  of  GASBUBL5  euid  produces 
an  output  tape  of  R,  R,  and  S array  multiplexed  as  in  GASBUBL5.  The 
data  are  equally  spaced  in  time  as  directed  by  control  cards.  The 
storage  foiroat  is  the  same  as  that  used  by  GASBUBL5  except  that 

(1)  only  three  arrays  are  multiplexed  (time  is  absent), 

(2)  ID(IOO)  is  3, 

(3)  XD(6o)  is  "GASRaiPL,"  and 

(4)  XD{42)  to  XD(59)  have  been  changed  to  reflect  the  nature  of 
the  data  stored  on  the  following  DR  (e.g.,  XD(57),  XD(58), 
and  XD(59)  are  all  set  to  the  time  step). 
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